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Abstract: We study the new boundary condition of the 0(n) model proposed by Ja- 
cobsen and Saleur using the matrix model. The spectrum of boundary operators and their 
conformal weights are obtained by solving the loop equations. Using the diagrammatic 
expansion of the matrix model as well as the loop equations, we make an explicit corre- 
spondence between the new boundary condition of the 0{n) model and the "alternating 
height" boundary conditions in RSOS model. 
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1. Introduction 

Boundary conformal field theories play an important role in many fields of theoretical 
physics, such as statistical mechanics, condensed matter or string theory. In order to study 
the properties of boundary conditions and boundary operators, it is useful to have at our 
disposal a microscopic description of the conformal field theories (CFT). The 0{n) model 
and solid-on-solid (SOS) models are the familiar examples which provide us with such 
a description of, in general irrational, CFTs with central charge c < 1. In the 0{n) 
model each lattice site is assigned an 0{n) spin, whereas in SOS models one associates 
an integer- valued height to each lattice point. In both theories, neighbouring sites are then 
coupled via suitable interactions. The heights are bounded from both sides in the so called 
restricted SOS (or RSOS) models; these models are known to describe rational CFTs with 
c < 1. 

Both the 0{n) and SOS models can be reformulated as loop gas models [1]. In this 
formulation, the 0(n) model makes sense for arbitrary real n and exhibit critical behaviour 
for \n\ < 2. The phase structure of these models is well understood. Interestingly, they are 
known to describe two CFTs of different central charges connected by a renormalization 
group flow. The loops behave differently in the UV (or dilute) phase and the IR (dense) 
phase. 

Some properties of the 0{n) and SOS models can be studied by putting them on 
a fluctuating lattice, i.e., coupling them to the two-dimensional gravity. The partition 
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function of such theories is given by summing up the partition functions of the model 
on all the different lattices weighted by their area. Actually, the 0{n) and SOS models 
on random - or dynamical - lattice are known to be described by the Feynman graph 
expansion of certain matrix models [1-7]. In this context, the continuum limit is achieved 
by tuning the potential couplings while sending the size of the matrices to infinity. In this 
limit one recovers the dynamics of the irrational CFT with c < 1 coupled to the Liouville 
gravity and reparametrization ghosts. 

The 0{n) and SOS matrix models are also useful in studying the conformally in- 
variant boundary conditions from the microscopic viewpoint. In this paper, we will be 
particularly interested in the boundary conditions of the 0{n) model recently proposed by 
Jacobsen and Saleur [8]. Instead of allowing the loops to touch the boundary freely, they 
weighted the loops touching the boundary differently from those which do not. They also 
considered the "L-leg" boundary operators on which L open lines end. The properties 
of such boundary conditions and operators were studied on a fixed annular lattice with 
L non-contractible loops introduced. They obtained a continuous spectrum of boundary 
operators and determined their conformal weights. The new boundary conditions were 
also put on a dynamical lattice by Rostov [9], where the correlation functions and the 
conformal weights of the L-leg operators were computed. The L-leg operators were also 
considered in some earlier works [10-12]. 

Another interesting fact is that the 0(n) model becomes equivalent to the RSOS model 
for some special values of n. In [8] it was proposed that the new boundary conditions of 
Jacobsen and Saleur correspond to the boundary conditions in RSOS model which force 
the boundary height to alternate between two values [13]. 

In this paper we study the property of boundary conditions and boundary operators 
of these models using the loop equations along the line of [9], but with more help of 
the matrix model formulation which is much simpler to handle than the combinatorics 
employed in the earlier work. We will be focusing on the dense phase, leaving the analysis 
of the dilute phase as a future work. 

The organization of this paper is as follows. In section ^ we introduce the new bound- 
ary conditions and boundary L-leg operators in the 0{n) matrix model following [9], and 
rederive the correlation functions and conformal weights of the L-leg operators from the 
loop equations. Then in section ^ we propose a description of the boundary conditions of 
alternating heights in RSOS matrix models. Using this we derive the spectrum of bound- 
ary operators as well as their conformal weights and correlators, again by solving the loop 
equations. Finally, in section ^ we show the equivalence of the 0{n) and RSOS models 
on discs by establishing a map between their Feynman graphs. We use this to derive some 
relations between the disc correlators, which are then shown to map the loop equations of 
one theory to the other. The last section |] is devoted to some conluding remarks. 

In Appendix |^ we review some basic facts on the Liouville theory approach to con- 
formal field theories coupled to two-dimensional gravity, and summarize the formulae for 
the conformal weight and gravitational dimension of the operators. Some detail of solving 
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the loop equation and reading off the gravitational dimension are given in Appendix ^ 



2. Boundary operators in the 0{n) model 
2.1 Definition of the model 

Let us consider a triangular lattice T with an 0{n) spin component associated to each site 
r, normalized so that tr Sa{r)St,{r') = SabSrr'- The partition function of the 0{n) model is 
defined by [14,15] 

Zr(T) = trn (l + ^5^^a(r)5,(r')j, (2.1) 

{rr') \ a / 

where T is called the temperature and the product runs over all links {rr') of T. Expanding 
the product into a sum of monomials, the partition function can be written as a sum over 
all configurations of self avoiding, mutually avoiding loops on F, 

Zr(T) = ^T"(''="s"')n#('°°P'). (2.2) 

loops 

Each loop is counted with a factor n, and the temperature T controls the average total 
length of the loops. When formulated in this way, the model makes sense for arbitrary 
real n. This model is known to exhibit a critical behaviour for \n\ < 2. Hereafter we 
parameterize n in terms of g or 6 as follows, 

n = — 2 cos (TT^f) = 2 cos (tt^^). (2.3) 

As a function of n, g is multi-valued. Different branches are known to correspond to 
different phases of the model [1]. 

The temperature controls the phase of the model. It is in the dilute phase at some 
critical temperature T = T*, and below T* it is in the dense phase. For generic n, the two 
phases are described by two irrational conformal field theories with central charges 

- - _ 

Cdense ^ -, ^ ; Cjiiute -L i i /i ' \^-^) 

1 — O i + U 

We will be focusing on the physics in the dense phase, which has the same behaviour as 
the fully packed loop model corresponding to T = 0. Hereafter we assume 9 E [0, 1] and 

= 1-9- 

Turning on the gravity corresponds to taking the sum over all the triangulated surfaces 
with a suitable weight, 

Z,y,iK,T) = ^/€-^(r)Zr(T). (2.5) 
r 

The parameter k controls the average area A(T) (the number of triangles) and is regarded 
as the bare cosmological constant. The continuum limit is obtained from the vicinity of the 
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critical line k, = k*{T) where the average area of the surface diverges. One can also allow 
the surfaces to have boundaries. For example, a disc partition function can be defined as 
the sum over the surfaces of disc topology, 

Zdy„(«,x,r) = -j^^^~^''^^x-''^^^Zr{T), (2.6) 

T: disc ^ ' 

where x is the boundary cosmological constant controlling the average boundary length 
(the number of edges along the boundary). In this case, we have to send x also to a 
critical value as k — > k*{T) so that the average boundary length diverges in the limit. The 
continuum limit is therefore parametrized by the renormalized couplings ji ^ k* — k and 
^ x—x*. Note that, when the disc has more than one boundary, a boundary cosmological 
constant may be introduced for each. 

Until recently, the only boundary condition studied in the 0{n) model was the Neu- 
mann boundary condition in which the spins at the boundary fluctuate freely. Based on 
earlier work [16, 17], Jacobsen and Saleur [8] proposed a new kind of boundary conditions 
in which the boundary spins are forced to take the first k of the n values. We call this the 
k-ih. JS boundary condition. Neumann and Dirichlet boundary conditions correspond to 
the special cases with k — n and k — 1, respectively. In the loop gas picture, the k-th JS 
boundary condition amounts to giving a weight k, instead to the usual n, to the loops that 
touch the boundary at lease once. Defined in this way, the JS boundaries make sense for 
non-integer k. 

Following [9], we consider the model on the disc with one Neumann and one JS 
boundaries connected by the boundary changing operators, 

Si= E Sa,---Sa„ §f= E Sa,---Sa,. (2.7) 

l<ai<---<ai<A; k<ai<---<a]^<n 

In the loop gas picture, they have L legs of open lines attached. They are called the blobbed 
and unblobbed L-leg operators [8]. They were named after the underlying Temperley-Lieb 
algebra though we will not need its detailed property in this paper. One of the important 
characteristics of these operators is that the lines from blobbed operators can touch the 
JS boundary whereas the line from unblobbed operators cannot. We will give the matrix 
model equivalent of these operators in the next subsection. 

2.2 The 0(n) matrix model 

The 0{n) matrix model is an integral over N x N hermitian matrices X and 1^, with a 
running from 1 to n. The partition function is given by [6] 



(2.8) 



JdxHdY^exp /3tr(-ix^ + ix^^-|j]i;2 + 5^Xi;^) 

a=l L 1=1 <J=1 

The Feynman graph expansion of Z generates all the dynamical lattices of arbitrary genus 
but without boundaries. The bare cosmological constant k, is given by 

(3^Nk^, (2.9) 
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and each loop formed by the propagators of Ya is multiplied by nT^^^™^^^\ Graphs of 
genus h are weighted by A^2-2/i^ ^^lat the planar graphs dominate the partition function 
in the large limit for a fixed k. Continuum limit is obtained by sending k k*{T) and 
iV — > oo in a correlated manner. 

The physics in the continuum limit depends on the temperature. Below the critical 
temperature T < T*, partition function is dominated by graphs with densely packed loops. 
Since the vertices with three legs of X do not play any role, one could study the dense 
phase using the definition ( |2.8| ) without the X^ term. Before proceeding, we make a slight 
redefinitions of the matrices X and Ya so as to simplify the integrand of ([2.8p, 



/n n 
dX^dYa exp -/?tr {v{X) XY^ 
a=l a=l 



(2.10) 



The potential V{X) is then given by 



y(x) = -^{x+\f-\[x+^)\ 



1 

3 



(2.11) 



As was mentioned above, generic quadratic potential V{X) could capture the physics in 
the dense phase. 

The disc partition function with Neumann boundary condition is given in the large N 
limit by 

$(a;) = -i(trlog(x-X)). (2.12) 

Because of the prefactor 1/(3, the leading contribution is independent of and the higher 
genus contributions are subleading at large A^. What will become more important later is 
its derivative, the resolvent 



(2.13) 



The derivative introduces one marked point along the boundary. One is supposed to take 
a; — > a;* = in the continuum limit [18]. The disc partition function with the k-\h JS 
boundary condition and one marked point is given by 



W{y) 



(2.14) 



where we suppress the fc-dependence of W for notational simplicity. To study the bound- 
ary changing operators, we also introduce disc two-point functions with one Neumann and 
one JS boundary conditions. 



Do{x,y) = ^ 



tr 



X 



(2.15) 
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We also consider the correlation functions of the L-leg operators, 



1 



Y" Y" 



(2.16) 



where the operators y| and are defined analogously to (P77|), 

{ai,--- ,ai}c{l,--- {ai,--- ,ai,}C{fc+l,--- 

The sums are taken over all different sets of L letters (so y\ consists of k\/L\ terms). 
These operators are the analogues in matrix model of the blobbed and unblobbed L-leg 
operators. 

In the following, we study the loop equations for the above correlators that are asso- 
ciated to Fa-derivatives. We will see that the Fa-derivative adds or removes one open line 
between boundary operators, so that the loop equations relate the correlators D^" with 
D'^f^. Among those correlators, d\ and Dq will be of particular importance because they 
can be determined from a closed system of shift equations. 

2.3 Loop equations 

We start from the loop equation which follows from the translation invariance of the mea- 
sure dYa- For any matrix F made of X and Ya, the following equality holds: 

If we introduce G = —{XF + FX), then F is formally expressed in terms of G as 

POO 

F= die^^Ge^^, (2.19) 
Jo 

Using them, the loop equation can be rewritten as 

Y \{w~. fdiie'^'Ge'"'),,) = (trGFa). (2.20) 

In the following we will apply this central identity to different G and derive some relations 
among our correlators d| and Dj^. 

2.3.1 Loop equations for Do and dI 

For later convenience, we begin by introducing the notation 

Hiy)= ^\ (2.21) 

y - Ea=i V 
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Let us first apply ( |2.2Q ) to G = e^'^HYa with a < k. Using the well known large 
factorization {trA tvB) ~ (trA) (tri?) and dropping the terms containing odd powers of 
in a correlator, we find 



Another relation can be obtained by applying ( [2 .201 ) to G = ^YaH: 

P{tre''''Y,HYa) = ^rff ((tre^^+^'^^F.i/y;) + (tre(^+^')^)) (tre^^i/). 

Now we make a Laplace transform with respect to i and £', using the relations 

[ d£e-^''tr(e^^ A) = tr( -^—A) , 
Jo \x-X J 

r e-^^'tr(e(^+^')^A)tr(e^^B) = trf -^a) * tr( -^b) . 
Jq \x — X J \x — X J 

Here we denoted by * the Laplace transform of the convolution 

dx' F(x') - Fix 



(2.22) 



(2.23) 



(2.24) 



F{x) * G{x) 



27ri X — x' 



G{-x') 



(2.25) 



and the contour of x' integration here encircles around the cut where F{x) has discontinu- 
ity. The loop equations ( |2.22[ ) and ( [2.23[ ) can then be rewritten into the form 



X 



(2.26) 



+ pWix)]*Doix,y). 



It only remains to sum over a from 1 to fc. Using the definitions ( |2.15| ) and (|2.16|) for Dq 
and dI as well as the equality 



k 

i J2 {tr{HY,'^^) ) = yDoix, y) - Wix) 

a=l 



(2.27) 



we obtain 



yDoix, y) - W{x) = Do{x, y) * (dJ (x, y) + kW{x) 
D\{x,y) = {kW{x) + D\{x,y)] * D^{x,y). 



(2.28) 



To obtain more useful equations for Dq and Dl, we subtract their non-critical parts 
and define 



do(x,y) = Do{x,y) - 1, 

kdi{x, y) = -Of (x, y) + kW{x) — y. 



(2.29) 
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It is natural to assume that do and di have the same cut in the a;-plane as that of W{x), since 
the cut is determined by the eigenvalue distribution of X and therefore does not depend on 
the correlators considered. The loop equation ( |2.28| ) can then be rewritten in terms of the 
discontinuity along the cut, 

do{—x,y)Discdi{x,y) +DiscW{x) = 0, 

1 (2.30) 
di{—x, ?/)Disc do{x, y) + ^Disc W{x) = 0, 

where Disc/(x) = f{x + iO) — f{x — iO). Using this one can show that the following 
quantity has no discontinuities and no poles in the complex x-plane, 

Pio(x, y) = d,{x, y)doi-x, y) + Wix) + jWi-x) - f (2.31) 

k k 

Moreover, its behaviour at large x is found from ( |2.29| ), 

do = -l + 0{-), d^ = -y + o{-), Pio = 0(-). (2.32) 
Therefore Piq should be identically zero. 



di{x,y)do{-x,y) + W{x) + jWi-x) -j = 0, 

k k 

do{x,y)d,{-x,y) + W{-x) + jW{x) - j = 0. 

k k 



(2.33) 



As compare to our loop equation ( |2.3(JD , the equation obtained in [9] (equations 3.13 
and 3.14) has one additional term coresponding to the JS boundary touching itself to break 
the disc into two pieces. Including such a term may be reasonable from the standpoint 
of the combinatorics because the JS boundary has fractal dimension 1/g. In the matrix 
model description, this boundary has classical dimension one and the term is missing. 
This discrepency comes from two different possibilities for defining the boundary Liou- 
ville potential. Because of the symmetry Ar^s(g) = Ag^ri^/g) between the dressed scaling 
dimensions, both point of view give the same scaling dimension for the boundary opera- 
tors. 

2.3.2 Loop equations for Dj^ and 

Let us next take G = e^'^YaH with a > k and apply ( |2.20D . Following the similar steps 
as in the previous subsubsection, we get 

Diix, y) = {n- k)W{x) * Do(x, y). (2.34) 

This is actually a special case of more general recursion relations between -Df_,_i and , 
and similarly between -D|^i and d\. To derive them, let us denote by {flj}, {hi] two 
arbitrary sets of order L + 1. They are both chosen to be subsets of {1, ■ ■ ■ , k} or {k + 
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1, • ■ ■ ,n} depending on whether we are interested in d| ot Dj^. Then we apply the loop 
equation ( [2 .201 ) to 

G = Y^,... _ \ ■ ■ • Y,,y, (2.35) 

where the derivative is with respect to ^a^+i, and sum over the sets [ai] and {hi}. The 
final result is 

Dl^,{x,y) = {k- L)Wix) * Dlix,y), 
Di+iix, y) = {n-k- L)W{x) * (x, y). 

These relations agree with the result of [9]. In terms of discontinuity along the cut, the 
loop equations become 

Disc Z^i+i(x, y) = {k- L)dI{-x, ?/)Disc W{x), 
DiscDj^^^{x,y) = {n - k - L)D-l{-x,y)D\^cW{x). 
The second equation can be extended to the case L = if one defines 

Di{x,y) = Do{x,y), Dl{x,y) = — (2.38) 

It was pointed out in [9] that these relations follow naturally if the Neumann and JS bound- 
aries are allowed to touch in dI but not allowed in Dq. 

2.4 Solution in the continuum limit 
2.4.1 The disc amplitude W 

To study the continuum limit, we introduce a small parameter e and set the unit lattice 
length to e. We define the renormalized bulk and boundary cosmological constants (/i, ^, Q 
by [1] 

e^;^ = K - K*, e^'^i = x-x\ eC = y- V* ■ (2-39) 

They blow up the neighbourhood of the critical point (k*, x*,y*) in the scaling limit e — > 0. 
Note that the JS boundaries have classical dimension 1 whereas the Neumann boundary 
has fractal dimension 1/g. The renormalized resolvent w(0 is defined by 

2V(.) uV'i-.) ^ (2.40) 

4 — 

In [1] the resolvent was obtained in the following parametric form, 

MS ^ 

^ = Mcoshr, = —— — cosh^fT. (2.41) 

Here M is related to the cosmological constant fi and the string susceptibility u via 



= 2gfi, u = ^^^^ = (2.42) 

As a function of ^, the resolvent has a cut along the interval ] — oo, — M] where the eigen- 
values of the matrix X are distributed. 
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2.4.2 The disc amplitudes and 

Here we wish to solve ( |2.33| ) in the continuum limit. To begin with, we need to find the 
critical value y* of the k-th JS boundary cosmological constant. Since it should not depend 
on X and k, one can determine it by requiring that cIq and di vanish at the critical point 

ix,y,K) = {0,y*,K*), 

y* = {k + l)W{0). (2.43) 
With a slight abuse of notations, we define the renormalized two-point functions by 

d,{x,y)=e-'/>^di{^,C), (2.44) 

where the scaling exponents ao, «i will be determined shortly. To the leading order in 
small e, the loop equations ( [2.33[ ) become 

di{^, C)rfo(-e, + ^(0 + 1^-0 - ^ = 0' 



(2.45) 



0. 



We dropped several terms in ( |2.33D such as polynomial terms in W{x) because they are 
subdominant for small e. By noticing rfj ~ for = ^ = and using the loop equations, 
one can determine the scaling exponents [9], 



ai 



1-9 -r9. 



(2.46) 



where r is related to k by 



k{r) 



sin (r + 1)tt9 



(2.47) 



sin rn9 

Hereafter we use r as the label of JS boundaries; it has a clear physical meaning as we will 
see later. We also express ( in terms of a new parameter a as 



CM 



M3 sm7i9 



cosh ga. 



2g sm7ir9 

Using (|2.41[), ([2.47[) and ([2.48D the loop equations can be rewritten as 



(2.48) 



di (r T y , ^) (r ± y , ^t) = CM^ cosh ^ ^ cosh ^ ^ , 



where 



C 



sin 7t9 



a 



{2r9 + 9-l). 



(2.49) 



(2.50) 



g sin 7r(r + 1)^^ 2 

We solve these shift relations in Appendix ^ using a slight generalization of [10]. The 
result is that do, di are given by the Liouville boundary two-point function [19]. See Ap- 
pendix for its explicit form. The scaling exponents of the correlators rfj oc are then 
read from their dependence on the boundary cosmological constant, and agree precisely 
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with ( |2.46D . This is enough to determine the gravitational dimensions of the boundary 
changing operators in the correlators do and di. 



= Ar,r, Af = A_,,_,_i, (2.51) 
where the gravitational dimension for the (r, s) operator is given by 

r-l-g{s-l) 



^r,s- ^ , (2.52) 

and is related to the conformal weight of the operators in CFT by KPZ relation ( |A.6p . See 
Appendix 1^ for more detail. 

2.4.3 The disc amplitudes Dj^ and dI 

Using the parameters (r, a) in the continuum limit and the equality 

w{t + ivr) — w{t — in) = sin tt9 sinh gr, (2.53) 

ig 



the loop equations ([2.37P can be rewritten into the form 



d\_^-^{t + ITT, a) — D|^-^(r — in, a) = — sin vr^ sinh (7ril'|(r, cr), 

(2.54) 



Dj^j^^ (r + ivr, cr) — -D/^+i (r — m, a) = sin nO sinh grDj^ v^'^)- 

'^9 



These shift relations take the same form as ( |A.12| ) satisfied by the Liouville boundary two- 
point function. So, in the continuum limit the correlators d\ and Dj^ are again given by 
Liouville boundary two-point functions up to factors independent of r, ex. The gravitational 
dimensions of the L-legs boundary operators are given by 

A| = A_,,_,_i, Af = Ar,r-L- (2.55) 



Through KPZ relation ( |A.6| ) this determines the conformal weight of the blobbed and 



unblobbed operators. The results are in complete agreement with [8] and confirm our 
identification of matrix model correlators with those of the 0{n) model coupled to gravity. 

As was noticed in [8], there is no reason for the parameter r to be quantized in the 
0{n) model. So we have a continuous spectrum of JS boundary conditions and the as- 
sociated boundary-changing operators. We conclude this section with a few remarks on 
some special values of r. First, for k = n ox r = 1 the JS boundary becomes the Neu- 
mann boundary. In this case the conformal weights of the boundary operators become 
= Si^L+i and 5^ = in agreement with the result [20] on flat lattice and [10] on 

dynamical lattice. Another interesting special case is the Dirichlet case, which corresponds 
to /c = 1 or 

r = (2.56) 
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In this limit, using the fact that the correlators of odd powers of Ya matrices vanish, one 
can prove the relations 



Dl{x,y) = yD^{x,y) - W{x) = W^' o; - X - V; / " 
The loop equation fovk — 1 then involves only one undetermined quantity, 

ydo{x, y)do{-x, y) + W{x) + W{-x) - y = 0. (2.58) 

We recovered the loop equation for the correlation functions of twist operators in loop gas 
model [10] up to normalization of correlators and parameters. Note that this simplification 
is a special feature of A; = 1 because the resolvents for the JS boundaries with A; > 1 
involve k non-commuting matrices Ya. 

3. Boundary operators in RSOS model 
3.1 Definition of tlie model 

In the RSOS height model [21], the local fluctuation variable (height) takes values in the 
integer set {1, • • • ,h — 1}. This model is also called A^-i-model, the integer set being 
identified with the nodes of the Dynkin graph for the A series. This graph is characterized 
by its adjacency matrix 

{1 if a ~ 6 , 
(3.1) 
otherwise. 

The indices a, b run over the nodes of the Dynkin graph, and a ~ 6 means that the nodes 
a, are linked. One can define the so called ADE-models in the same way from the Dynkin 

graphs of the ADE Lie algebras [22] . 

The RSOS model on a fixed triangular lattice with possible curvature defects is de- 
fined as the statistical sum over all the height configurations. Each height configuration 
is weighted according to the following rule [23]. To each site of height a one assigns the 
local Boltzmann weight 

Wo{a) = Sa, (3.2) 
and to each triangle with the heights a, 6, c at the three vertices one assigns 

WAia, b, c) = -j=5abhc + ^ i^abGbc + hcGca + ^caGab) ■ (3.3) 



Sa T\/Sa 

Here T is the temperature, and -S'a are the components of the Perron-Frobenius vector 

2 . /7ra 
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which is the eigenvector of the adjacency matrix with the largest eigenvalue 2 cos (7r//i). 
The weight ( |3.3| ) in particular requires that the heights of any two adjacent sites can differ 
at most by a unit. Thanks to this, the height configurations can also be described by the 
contour lines (loops) along the edges of the dual lattice [15]. The average total length of 
the loop is then controlled by the temperature. 

The phase diagram of the RSOS model is the same as that for the 0{n) model; it 
is in the dilute phase at some critical temperature T*, and in the dense phase at lower 
temperatures. Since the lattice is filled by the loops in the dense phase, one can study this 
phase using the Boltzmann weight ( |3.3D without the first term. 

Conformal boundary conditions are realized microscopically as suitable restrictions 
on the heights on the boundary. In the dense phase there are two kinds of boundary condi- 
tions. The "fixed" boundary condition a requires the boundary to take the constant height 
a. The "alternating" boundary conditions (ab) require the sites at the boundary to take the 
two adjacent heights a and b alternately, like ababab...ab. Bauer and Saleur [13] studied 
both types of boundary conditions on the flat lattice, and found that conformal weights of 
the boundary changing operators ^i?" and ^i2>^(ab> ^j-g given by 5i^a and 6a,i- In the follow- 
ing we will show that the boundary operators of conformal weight Sr^s can be realized as 
sj^(rr+i) vv^hen r > s and when r < s. 

3.2 The RSOS matrix model 

The RSOS matrix models [24] are the simplest examples of the AD E -matrix models. 
The fluctuating variables are /i — 1 hermitian matrices Xa associated to the nodes a = 
1, • • • ,h — 1, and h — 2 rectangular complex matrices Cab = CL associated with the 
oriented links (ab). Xa has the size Na x A^^^ while Cab has the size NaX Nh. The partition 
function is given by the integral 



J] dXa n dCab exp ( - (3S[X, C]), 

a (ab) 

S[X, C]=J2 Satr (l-Xl - \xA + ^ tr (^CabC,a - XaCabQ 

^ ^ (ab) 

where {ab) runs over oriented links of the Dynkin graph. 

We take the limit of large P, large Na keeping their ratio fixed, 

^ = ^' = fixed. (3.6) 
The constant k plays the role of the bare cosmological constant. Perturbative expansion 



of Z gives the sum over height configurations with Boltzmann weights ( |3.2| ) and (p^, but 
now on dynamical lattice. The planar graphs dominate the partition function in the large 
P limit, and the higher genus terms are suppressed by powers of 

One can view the system as the gas of loops on dynamical lattice which are formed by 
the propagators of Cab, Cba and separating the domains of heights a and b. The temperature 
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T is regarded as the fugacity for the total loop length. Again, the term in S[X, C] can 
be dropped when one is interested in the dense phase. 

Also, the Feynman rule is such that each connected domain of height a gives rise to a 
factor 

(3.7) 

where x is the Euler number of that domain. For disc graphs, this rule amounts to assigning 
the factor Sa to each "outermost" domain of height a touching the boundary, and Sa/ Sb to 
each domain of height a surrounded by that of height b. 

Before proceeding, we make a certain linear change of matrix variables to rewrite the 
partition function as 

Z = jUdXall dCab exp ( - /5 5^ SatrViXa) +(3^2 XaCabCba) , (3.8) 

a (ab) a (ab) 

where the potential V{x) is the same as the one ( |2.1 ID for the 0{n) model. As was men- 
tioned above, generic quadratic V can describe the physics in the dense phase. 
We will consider the following correlators 

'^aix) = -^(trlog(x-X„)), 

^ 1 (3.9) 

'^{ab)iy) = -^{tr\0g{y - CabCba))- 

They correspond respectively to the disc partition function with the boundary conditions a 
or (ab). The more important quantities in the following analysis are their derivatives, the 
loop amplitudes 

d IX 

Wa(x) = - — <i>a(x) = -(tr 



The loop amplitudes Wa{x) axe the Laplace images of the loop amplitudes with fixed 
boundary length, 

iy,(£) = -^(tre^^»>. (3.11) 

To study the spectrum of boundary operators, we also consider the disc two-point functions 
with a and (be) boundary segments, 

= 1(„(_J_S^,_2^SL)). (3.12) 

Here 

^ab = CadiCdid2 ' ' ' ^d^^ib (3.13) 
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Da,Mii,y) = -(trfe^^^S^,— ^^SL)), (3.14) 



is the product of C"s along the shortest path from a to b, so that each height di appears 
only once. We also assume that the path does not contain the node c. There are L open 
contour lines stretching between the two boundary changing operators S^^. Their Laplace 
transform, 

will also appear in the derivation of the loop equation. 
3.3 Loop equations 

Here we derive some loop equations for disc correlators by using the same technique that 
was used in the 0{n) model. Let us start from 

= -(trF(aaX, + X,aa)>, (3.15) 

P \ OUab ij I 

where F is a A^^ x A/";, matrix made of X and C matrices. By inserting 

die^^^Ge^^\ (3.16) 

^0 

we obtain the following equation 
3.3.1 Recursion relation for £'a,(fec> 

To begin with, we derive a loop equation involving Da^i^bc) assuming a 7^ 6 and that c is not 
on the shortest path connecting a and h. Let d be the node adjacent to a along that shortest 
path, so that S^f, = Cad^db^ and = Sj^^^Cda- Applying ( FTVD to 

G = e^'^-a. S^,-^ ^ S^-^ , (3.18) 

gives the following recursion relation 



Da,ibc){i',y)= / diWa{i + nD,^^bc){i,y). (3.19) 
Its Laplace image reads 

Da,(bc){x,y) = Waix) * Dd,{bc)ix,y), (3.20) 

with the * product defined in (|2.25D . In this derivation, it is important that Cad appears 
only once. The relation ( [3.20D can also be obtained by an explicit integration over the 
Cad matrix using the Gaussian measure. From the viewpoint of height configuration, this 
can be understood as cutting the disc into two pieces along the contour line that stretches 
between the two boundary operators and separates the domains of height a and d. See 
Figure [1]. This recursion relation allows us to express all the disc correlators Da,{bc) in 
terms of the 2h — 2 basic correlators -Da,(aa±i)- 
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Figure 1: The disc graph contributing to Da^[bc) can be decomposed into two discs along the 
contour line separating the domains of height a and d and connecting the two boundary operators. 
^a,{bc) is therefore written as a *-product of Wa and -Dd,{bc)- 

3.3.2 Bilinear functional equation for Da^(^a.b) 

Let us next apply ( |3.17D to 



^ = 6^'^" ^—^Cab. (3.21) 

Following the similar steps as before we obtain the loop equation 

POO 

y / rf££),,(,fe)(£ + f,y)A,(,„)(£,y)=y£),,(,,)(f,|/)-W^,(0. (3.22) 





After the Laplace transform with respect to i, it becomes 

Da,{ab){x,y) = ^Wa{x) + Da,{ab){x, y) * Db,(^ba) {x , y) . (3.23) 

This relation can also be understood from the viewpoint of the height configurations. Re- 
call first of all that the {ab) boundary emanates contour lines separating the heights a and 
b. Let us then think of a graph participating in Da,(ab), and cut it along the contour line 
whose endpoint is the closest to the right end of the (ab) boundary. There may be no such 
contour line because the (ab) boundary may have zero length; in such a case the graph 
has constant boundary height a and contributes to Wa{x). Otherwise the graph can be 
decomposed into two pieces contributing respectively to Da^i^ab) and Db^{ba)- See Figure ||. 
These two possibilities correspond to the two terms in the RHS of ( |3.23| ). 

The disc correlators Da^i^ab) and Di,^(ba) satisfy two relations, namely ( p.23| ) and another 



equation obtained by exchanging a and b. To get the more useful equation we take the 
discontinuity of these equations along the cut in the x-plane, 

DiscDa,{ab){x,y) = -DiscWa{x) + Db,{ba){-x,y)DiscDa,(ab){x,y), 

\ (3.24) 
T>hcDb^(ba){x,y) = -DiscWb{x) + Da,{ab){-x,y)DiscDb,(ba){x,y). 

Using these and taking into account the large-x asymptotics of Wa and Da,{ab), we find 
yDa,{ab) {x, y)Db^{ba){-x, y) 

-yDa^t^ab){x,y)-yDb^(ba){-x,y) + Wa{x)+Wb{-x) =0. (3.25) 
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y 



[ab) 




Figure 2: There are two types of graphs contributing to -D^ (^6) • The graph may have (ab) boundary 
of length zero so that the whole boundary has height a. If {ab) boundary has nonzero length, one 
can cut it into two pieces contributing to D^^i^ab) ^rid -Db,(6a)' respectively. 



By introducing 

da,{ab){x,y) = ^{Da,(^ab)ix,y) - l) 

we can rewrite the equation ( p.25D into a simpler form, 

da,(ab){x,y)db,{ba)i-x,y) + Waix) + Wb{-x) = y. 



(3.26) 



(3.27) 



Here the normalization of d^^afc) is somewhat arbitrary, and we chose the symmetric nor- 
malization for simplicity. 

3.4 Solution in the continuum limit 

To study the continuum limit, we introduce a small parameter e and the renormalized 
couplings (/i, ^, C) in the same way as we did for the 0{n) model, 

e^ix = K- k\ t^/^i = x-x\ < = y- y*. (3.28) 

In [1,7] the loop equations for Wa{x) has been solved under the natural ansatz, 

Wa{x) = SaW{x). (3.29) 



W{x) was then shown to satisfy the loop equation ( [2.1 3D for the resolvent of the 0{n) 
matrix model with n = 2cos(7r//i). So we borrow the solution from the 0(n) model 
under the identification 

e=l-g = ^. (3.30) 



We define the renormalized resolvent w{^) in the same way as in ( [2.40[ ). The solution in 
the continuum limit is given by 



i = M cosh r, 
Again, M is related to ji via (|2.42|). 



M3 

= -— cosh gr. 



(3.31) 
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Now let us solve the loop equation ( [3.27D for the correlators da^i^ab)- First we need to 

determine the critical value of the (a6) -boundary cosmological constant y* . We require 
that da^i^ab) and vanish at the critical point (/t*, x*, y*) and find 

y* = Wa{^) + W,{<d) = {Sa + S,)W{<d). (3.32) 

Next we renormalize the disc correlators near the critical point as 

1 (3.33) 

The loop equation in the limit of small e is 

da,{ab){^,Odb,{ba){-^,0 + SaWiO + SbW{-0 = C- (3-34) 



By substituting ( 3.31 ) and 

Q = cosh ga (3.35) 

2^ 



into ( [3.34D , the loop equation finally becomes 

da,{ab) {r - in/2, (T)(ib,(fea) {r + i7r/2, a) 

= ^'^^ cosh (^ 3{r + a) + aab ^ ^^^^^ j^^(r--aH^ j (3.36) 



9 

where Uab for 6 = a ± 1 is given by 



m f a + 6\ 
Uab = ±— 1 ; — = -aba- (3.37) 



2 V h 

The above equation has the same structure as ( |2.49| ), so the solution is given in terms 
of the Liouville boundary two-point functions. We can now use the formulae in Appendix 
^ and and determine the gravitational dimensions of the boundary changing operators 

Aa^^ab) = Aa^a, Afc,(fc„) = Ai^b-b whcn 6 = 0+1, 

Aa,{ab) = Ai_a-a, Ab_(6„) = Afc^fo whcn b = tt - 1. (3.38) 

This implies that the boundary operators and '^B^"'^'^^ in RSOS model have con- 

formal weights Sa,a and 5a-i,a5 in full agreement with the result of Saleur and Bauer. 

From the comparison of the loop equation ( |3.20D with ( |A.12D , it follows that the cor- 
relators c?a (be) are all given by Liouville boundary two-point functions. The operators 

a-Lj^{aa+l) a+Lj^{aa-l) j-j^gj^ shown tO haVC COnformal weights 5a,a+L and 5a-l,a+L, 

respectively. By varying the integer parameters a and L, the whole spectrum of boundary 
operators for this rational CFT is recovered. This result proves the conjecture of [13] on 
the scaling dimensions of boundary changing operators in the RSOS model. 
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4. The map between the two models 



It has been known for a long time that the 0{n) model and SOS models are described by 
the same class of conformal field theories. In particular, the RSOS model with h — 1 nodes 
was known to have the same partition function on the plane as the rational 0{n) model 
with n = 2cos(7r//i). Here we wish to explore this correspondence further, focusing 
mainly on the theories on the disc. 

As was used in the previous section, the 0(n) model on the disc with Neumann bound- 
ary condition is equivalent to the RSOS model with fixed-height boundary condition. The 
disc partition functions are related via 

Waix^SOs) = SaWiXoin))- (4.1) 

In [8] it was shown that the (a6)-type boundaries of RSOS model correspond to the JS 
boundaries of the 0{n) model labeled by integer r. The annulus partition functions of 
the 0{n) model with Neumann- JS boundary conditions were shown to agree with those 
of RSOS model with a- {be) boundary conditions. It was also noticed there that one needs 
to introduce L non-contractible loops on the annulus of the 0{n) model, corresponding to 
the distance between a and {be) labelling the two boundaries of the RSOS model. 

Following their idea, we wish to relate the disc correlators of the two models on dy- 
namical lattice. More explicitly, we will find out the relation between d|, Dj^ of the 0{n) 
model and Da+L,{aa-i), Da-L,{aa+i) of the RSOS model. Our derivation of the relation 
is based on the loop equations and the diagrammatics of the matrix models, and does not 
rely on the continuum limit. 



The first step is to relate the boundary cosmological constants. From the relation ( |4.1D 
we simply relate the x's by 

a^RSOs = a;o(n) (4.2) 

for Neumann boundary of the 0{n) model and the fixed height boundary in RSOS model. 
Similarly, it follows from ([2.43|) and (p.32p that the y's for the r-th JS boundary and the 



{ab) boundary should be related via 



yo(n) = {{a,b) = (r,r + l) or {h - r, h - r - 1)) , 

yo{n) = {{a,b) = (r + l,r) or {h - r - l,h-r)). 

Oh 



(4.3) 



4.1 Relations between Feynman graphs 

Instead of finding the correspondence of disc correlators quickly from loop equations, 
let us explain how the correlators of the two matrix models should be related from the 
viewpoint of Feynman graph expansion. 

The underlying idea is very simple. Each Feynman graph of the RSOS matrix model 
describes a dynamical lattice with a height assigned to every face, so that one can draw 
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contour lines separating the domains of adjacent heights. If we focus only on the contour 
lines and forget about the heights, then what we get is nothing but the Feynman graph of 
the 0{n) matrix model. We thus compare each Feynman graph of the 0{n) matrix model 
with the sum over all the Feynman graphs of the RSOS model having the same contour 
line configuration but different height assignments. 

4.1.1 Resolvents 

To begin with, let us consider the relation between the resolvents of the RSOS and the 
0(n) matrix models, 

Wa{x) = SaW{x). (4.4) 

Each graph contributing to Wa{x) has the unique outermost domain of height a, and is 
inscribed by several subdomains of height a ± 1 . Each subdomain may be inscribed by 
several subdomains of adjacent heights, and by iterating this a finite number of times one 
can cover the entire disc. Now let us sum over all the height assignments in the interior. 



Using the rule explained after (13. 7D, we first assign Sa to the outermost domain of height 
a. To each of its subdomains one can assign the height a + 1 or a — 1, which gives rise to 
a factor 

Sa+l , Sa-1 o /-A i:\ 

+ = 2 COS - = n. (4.5) 

By repeating this and going step by step to the interior, one ends up with the Feynman 
graph of the 0{n) matrix model with a factor n assigned to each loop. This explains the 
relation ( ]4.4D at the level of the Feynman graph sum. 

4.1.2 Disc correlators and £'a,(afe> 

Next, we use a similar argument to show the relation 

Da,{ab) (x, y) = Dq (x, y/ Sa) , (4.6) 

where we assume b = a + 1 for simplicity, and the JS boundary is expected to be labelled 
by r = a from (|4J|). To show this, we consider Feynman graphs contributing to the LHS 
with the power l/y™"^^. Such graphs have the (a6) -boundary of length m, and there are 
therefore m open contour lines ending on the (a6) -boundary. If we cut the graph along 
these contours, it would decompose into rua pieces of boundary height a and mi, pieces of 
boundary height b, where rUa + nib = m + 1 and ma > 0. This means that the graphs have 
m + 1 outermost domains of heights a or b. 

We now take the sum of such graphs over different height assignments in the interior 
but for a fixed contour line configuration, to obtain a graph of the 0{n) matrix model. Each 
contour line in the interior is assigned a factor n in the same way as before. Collecting the 
factors associated to the outermost domains we find, 

^-rn-l^rn^^rn, = (S./y)"^^' X (4.7) 
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Figure 3: The graphs contributing to -Da,(afe) have open contour lines separating the domains of 
height a and b. By shrinking the b part of the (a6) -boundary, those contour Unes turn to form loops 
touching the boundary. 

We then deform the contour line configuration to a loop configuration by shrinking 
the 6-part of the (a6) -boundary as shown in Figure The m open contour lines then 
turn into rrih closed loops touching the boundary m times in total. The expression for the 
weight ( |4.7p then implies that the (a6) -boundary of the RSOS matrix model is mapped to 
the k-th JS boundary of the 0{n) model, with boundary cosmological constant y/Sa and 

. Sa+i sin7r(a + 
= = ■■ 77 — ■ (4.8) 



Thus we have shown (P3|). It also implies a = r in agreement with (|4.3|). 



4.1.3 Disc correlators Df and (ba) 

Using the same argument, let us next show the equation 

Si 

y 



Db,(ba){x,y) = ^(W{x) + ^Dl{x,y/Sa)). (4.9) 



We again focus on the graphs contributing to the LHS with power 1/?/'"+^, which have m 
open contour lines ending on the (6a) -boundary. Such graphs have rua outermost domains 
of height a and outermost domains of height b, with the condition 

ma + mh = m + l, > 0. (4.10) 

We perform the sum over the height assignments in the interior and map the graphs to 
those of the 0{n) matrix model. 

By shrinking the 6-part of the (6a) -boundary, we get the graph in the 0{n) matrix 
model which generically has one open line connecting the two boundary changing opera- 
tors in addition to m?, — 1 loops. They altogether touch the JS boundary m — 1 times in 
total. It is important to notice that the open line can end on the JS boundary. We should 
therefore identify the boundary operators as the one-leg blobbed operators . The situ- 
ation is illustrated in Figure 0. The graph of the 0{n) model we thus obtained has the 
following weight from the outermost domains. 
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Figure 4: The graphs contributing to -D;, (^a) have open contour lines separating the domains of 
height a and b. By shrinking the b part of the (6a) -boundary, those contour lines turn into some 
loops touching the boundary and a line connecting the two boundary-changing operators. 



The same graph and weight can be obtained from the Feynman graph expansion of the 
second term in the RHS of (|4-.9D. Note that the additional factor is inserted because 
the line connecting the two Yf can be made from propagators of Yi, ■ ■ ■ or Yk, leading to 
a factor k. The first term of the RHS, on the other hand, is the sum over the exceptional 
graphs of the RSOS model corresponding to m = 0, namely, those graphs which have the 
(ba) boundary of length zero. The sum over such graphs is simply the leading term in the 
1 / ^/-expansion of the LHS and therefore given by Wb{x) /y. This finishes the diagrammatic 
proof of 

The relations ( |4.6| ) and ( p!9| ) can be used to show that the loop equations of the 0{n) 
matrix model ( [2.33[ ) and the RSOS model ( |3.25| ) are mapped to each other. These relations 
also explain that the operators '^^(f'a+i) and "+i5(a+ia> have the same conformal weight 
as that of Sq and sf between the Neumann and the A;(a)-th JS boundaries. 

4.1.4 Disc correlators of L-leg operators 

It is obvious how to extend the correspondence to the disc correlators of L-leg operators 
using the argument of summing over height configurations and shrinking the 6-part of the 
(ab) boundary. We skip the details and present the final results. 

{n-k-L)\ 



Da-L,(aa+l){x,y) 
Da+L,{a+la){x,y) 



Sa-1 ' ' ' Sa-L 
Sa+1 ■ ■ ■ Sa+L 

y 



(n 



k)\ 



-Dtix,y/Sa) 



Wl(x) 



{k-L)\ 



k\ 



Dlix,y/Sa] 



(4.12) 



where Wl{x) is determined by the recursion relation 

DiscWL+iix) = WLi-x)DiscW{x), 



Wi{x) = W{x) 



(4.13) 



These include the results of the previous subsubsection as special cases. It is also easy 
to show that they relate the recursion relation of the 0{n) model ( |2.36D to that of the 
RSOS model ( |3.2Q ). They also explain that the operator C'-Lj^iaa.+i) j^^g ^j^g same confor- 
mal weight as §f between the Neumann and the /c(a)-th JS boundaries, and similarly for 

a+L^{a+la) gll_ 

The first term in the second line of ( |4.12| ) is equal to the leading term in the 1/y- 
expansion of the LHS, and the corresponding graphs have {a + la) boundary of zero 
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length. Therefore, Wl{x) is a disc one-point function of a boundary operator with L 
nested loops attached, corresponding to the fusion product of two L-leg operators. Such 
an operator should be described in terms of star operators [11]; the star operator 5^ is a 
source of L open lines and is allowed to exist between two Neumann boundaries. In [1 1] 
its gravitational dimension was found to be Ai^^+i, and the disc two-point functions were 
computed in the continuum limit. Our Wl can be calculated in the same way. Using the 
standard parameterization ^ — M cosh r in the continuum limit we find, 

L . L 

■j^J-^ smTTg[k + 1) ^ ^ 

5. Concluding remarks 

In this paper we studied the new boundary condition of the 0{ri) model proposed by 
Jacobsen and Saleur. By using the matrix model formulation, we were able to relate them 
to the boundary conditions of RSOS model with alternating heights. The loop equations 
turned out to be a very efficient tool in calculating the spectrum and the conformal weights 
of boundary changing operators. 

Our techniques based on matrix model and loop equations are applicable to the anal- 
ysis of more involved situations, such as discs with several JS boundaries labelled by dif- 
ferent k. An interesting problem is to study the spectrum of boundary operators between 
two JS boundaries. (On a regular lattice this has been done in the recent work [25].) 

Another natural and interesting question will be to ask how our results can be extended 
to the dilute phase. Since the lattice will no longer be packed densely by the loops, one 
would expect a conformal boundary condition for which some sites on the boundary have 
no open line attached. We will therefore need to generalize the JS boundary so that it can 
have vacancies. It will be an interesting problem to study the renormalization group flow 
for the fugacity associated to the vacancy. We hope to address this issue in the future. 
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A. Gravitational dressing, Liouville theory and KPZ 



Here we summarize some basic facts on conformal field theories coupled to gravity and 
boundary Liouville field theory. More details can be found in [11, 12] and the original 
paper [19,26,27]. 

Let us consider the 'matter CFT' of the central charge 

c=l-^<l. (A.1) 

When 9 = 1/h for a positive integer h > 3, the theory is rational and corresponds to 
the minimal model of the unitary series {h — l,h). Turning on the gravity corresponds to 
summing over different metrics and topology of the two-dimensional space. After gauge 
fixing it amounts to coupling the CFT to the Liouville field and be ghost system. 

Let us take the matter CFT to be the Coulomb gas model [14, 15] described by a scalar 
X- The vanishing of the total central charge puts a certain condition on the matter and the 
Liouville background charges eo and Q. We can parameterize them as 

g = ^ + b, eo = ^-b. (A.2) 

Here b = y/l — 9 is the Liouville coupling constant. 

In this model we consider the matter field e*^'^ "^ of conformal weight 6r,s- With a 
suitable gravitational dressing it becomes an operator of conformal weight one, 

T(2bPrs) 

Br,s = — exp {ier,sX + (^r,s(t>) ■ (A.3) 

TT 

Here various parameters are related as follows, 

and (r, s) is a pair of positive integers labelling degenerate representations of the matter 
CFT. The matter conformal weight 5^,5 is given by 

(r/b-.b)^-(l/b-b)^ 

Or,s = ^ ■ (A.5) 

Introducing gravitational dimension A = (2P — eo)/2b and string susceptibility 7str = 
—9/{l — 9), one can write the KPZ scaling formula [28-30], 

^^A(A-^ (A.6) 

1 - Tstr 

Note that there is another way of gravitational dressing, /9 = (5/2+|P|,as was considered 
in [9]. As an example, the boundary identity operator can be dressed by e'^/^ instead of 
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gfa0 jj^jg suggests that there are two boundary cosmological couplings, and one has a 
fractional dimension with respect to the other. 

After turning on the gravity, correlators no longer depend on the positions of the op- 
erators inserted because one has to integrate over the positions of those operators. The 
dimensions of the operators therefore cannot be read from the position-dependence of 
their correlators. Instead, they should be read off from the dependence of correlators on 
the cosmological constant /i. If we restrict to discs, then the amplitudes with n boundary 
operators Bp^ and m bulk operators Vkj scale with ji as 

{^^Bpf- ■ ■ ■ ^-Bpi^ Vk,--- Vk^) oc fi\ (A.7) 

with 

7 = (1 - " - - + l^'l + I'^'O- 

i = l j = l 

Of course, the gravitational dimensions of the operators can be read more explicitly from 
the more detailed form of the amplitudes. 

As an example, let us consider a disc with two boundary segments, labelled by bound- 
ary cosmological constants (^i and ^2 and connected by the operators '^'^B'^^ . In boundary 
Liouville theory, boundary cosmological constant ( is the coefficient of the boundary in- 
teraction e*"^. Following [19] we use a parametrization of ( similar to ( |2.48p . 



C= a/— ^coshbV. (A.9) 

The computation of the disc amplitude involves the disc two-point functions of the Liou- 
ville and matter CFTs. The Liouville and matter part of the correlator factorize, and the 
matter part gives only a C-independent factor. The Liouville part is given by [19] 

i^^Bp^^ix) ^^Sp^KO)>„e = A{P)d{\Pln, r2)|x|-2^«3-/3) ^ (A.IO) 

with 

, . dijj / cos cjTi 008 07X2 sinh27rPtt;/b 2Pb\ ^ 

Ind P,ri,r2 = - / — — — — . (A.ll) 

J _^ uo \ smh TTCJ smn vrct;/ txuj J 

Here A{P) is a known function of P and is related to the "leg factor" arising from different 
normalization of the wave functions. It is independent of r's and therefore unimportant. 
On the other hand, the r-dependent part ( |A. 1 1| ) is expressed in terms of the double-sine 
function of pseudo-periods b and 1/b [31]. It satisfies an important shift relation involving 
both P and r^, 

d{P,Ti + m,T2) - d{P,Ti-iTi,T2) oc sinh bVi (i(P - b/2,ri,r2), (A.12) 

up to a r-independent factor. Shifting P by — b/2 corresponds to changing the label of the 
operator from (r, s) to (r, s + 1). 
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B. Solving the loop equation 

Here we solve the loop equation ( |2.49| ) 

rfi (r T y , a) rfo (r ± y , a) = CM^ cosh ^ cosh ^-^ ^ 

We define the function ^^(r, cr) by 



c^a(r,a) = ^M'^«expM,(r,a) (B.l) 

with + = (yf, and denote by Uaiy^-, c) their Fourier transform with respect to r. 
Let us take the log and the Fourier transform of the loop equation. Using 

dre'-^ log (2 cosh ) = . , , , (B.2) 

the loop equation for Uai^^, o") becomes algebraic, 

e^^..i(c.,a) + e^^^o(^,^) = -^_^^^e±-/^. (B.3) 

CO smh^vrcij/gfj 

Solving this in favor of iia and Fourier transforming back, we find 

, , f°° du cos auj cos TLJ , ^ /nuj , iauj\ 

Ua{T,a) = — — -smh — ± , (B.4) 

uo smh TTo; smh TTu;/ (7 V 2 g / 

where plus sign is for ui and minus sign for uq. One recognizes the same functional form 
as the Liouville boundary two-point function ( |A.l 1[ ). 

By comparing ( |B.4| ) with ( |A.l 1[ ) one finds the value of P for the boundary-changing 
operators. Then by using the scaling law ( |A.8| ) one can determine the exponents and ai 

2bPo = ao = re, 2bPi = ai = l-e-r9. (B.5) 

Another way to find is to analyze the two-point function at C = = 0, 

ci(^,( = 0,/i = 0) ~e"»l"l. (B.6) 

Setting a = iix /2g and r ^ 00 in ( |B.4| ), the dominant contribution to the integral is from 
the vicinity of the second order pole at = 0. Using 

due'^^ , , 

^ = -7rr, (B.7) 



one finds Ua — 2b P\t\ and recovers ( |B.5[ ) again. 
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